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Abstract. Two-particle correlation functions are calculated for bosons emitted 
from a localized thermal source (the "glow" of a "hot spot"). In contrast 
to existing work, non-equilibrium effects up to first order in gradients of the 
particle distribution function are taken into account. The spectral width of 
the bosons is shown to be an important quantity: If it is too small, they do 
not equilibrate locally and therefore strongly increase the measured correlation 
radius. 

In memoriam of Eugene Wigner and Hiroomi Umezawa. 



1. Introduction 



In the field of relativistic heavy-ion collisions the analysis of Bose-Einstein corre- 
lations j|] has attracted much attention recently. The general hope is to extract 
information about the size of a source radiating mesons by studying their two- 
particle correlation function ^. These correlations are typical quantum effects, 
hence quantum field theory is a proper framework to describe the problem the- 
oretically. Such efforts have been undertaken with great success for many years 
[H H; Hj 0|j with a big emphasis on the quantum properties of mesons that were 
emitted from a distribution (in space and time) of classical currents. The merger of 
these theoretical considerations with space-time distributions of quasi-particles gen- 
erated by simulation codes for relativistic heavy ion collisions has led to predictions 
of correlation radii that are in rough agreement with experimental data. 

The problem of a free quantum field radiating from a classical current is ex- 
actly solvable pp.438], also at finite temperature |9|. However, in the original 
field where Bose-Einstein correlations have been used to measure source radii, i.e., 
in the Hanbury-Brown-Twiss analysis of star" light" , as well as in the analysis of 
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relativistic heavy-ion collisions, the radiation source is a localized thermal distri- 
bution. Such a state is a non-equilibrium state, hence somewhat difficult to handle 
theoretically. 

As was pointed out above, the usual way to circumvent this problem is the ap- 
proximation of a non-equilibrium state as the superposition of an infinite number 
of classical currents, with a certain current distribution in space and time (mostly 
assumed to be gaussian) . From the viewpoint of quantum field theory this ap- 
proximation is somewhat unsatisfactory, if not questionable: True non- equilibrium 
effects on the propagating particles are neglected, only local equilibrium effects are 
maintained. 

Since also in relativistic heavy-ion collisions the radiating system is off equi- 
librium, one question arises immediately: Are there non-equilibrium effects on the 
quantum Bose-Einstein correlation ? A secondary question is, whether the corre- 
lation effects generated by a distribution of classical currents and those generated 
by a thermal source are equivalent or not. 

To address these questions we consider a stationary temperature distribution 
with limited spatial extension, and calculate the two-particle correlation function 
of bosons emitted thermally from such a distribution. In other words, we are 
investigating correlations in the "glow" from a hot spot - in close analogy to the 
physical situation present in astrophysics, but on the quantum level. 

To this end we formulate the two-particle correlation problem in a field theoret- 
ical method suited to handle non-equilibrium states. However, since the solution of 
the full problem is beyond our capabilities, we restrict ourselves to the same order 
of accuracy that is reached in standard transport theory: Our correlation function 
incorporates non-equilibrium effects beyond the local equilibrium function, but only 
up to first order in the gradients of the "temperature" distribution of this system. 
Furthermore, the discussion is limited to a static non-equilibrium system. 

In another sense we treat the correlation problem more consistently than in 
standard transport theory: We take into account a nonzero spectral width for the 
particles we consider. This is necessary, because at nonzero temperature every 
excitation acquires a finite lifetime due to collisions with the medium We 
describe this finite lifetime by attributing a certain spectral width 7 > inside the 
hot spot also to asymptotically stable particles. For strongly interacting particles, 
like e.g. pions, we may assume that such a spectral width is due to the coupling to 
A33-resonances. 

Quantum field theory for non-equilibrium states comes in two flavors: The 
Schwinger-Keldysh method |ll[ and thermo field dynamics (TFD, see |l^ for 
references to the original work). For the purpose of the present paper, we prefer 
the latter method: The problem of an inhomogeneous temperature distribution has 
been solved explicitly in TFD up to first order in the temperature gradients 
This solution includes a nontrivial spectral function of the quantum field under 
consideration. It employs a perturbative expansion in terms of generalized free 



fields with continuous mass spectrum pO[ |. 

The paper is organized as follows. The next section contains a brief introduction 
to the formalism of thermo field dynamics for spatially inhomogeneous systems. 
In section 3 we derive expressions for the two-boson correlation function in non- 
equilibrium states as well as its local equilibrium approximation. Section 4 contains 
a study of the effect that is exerted on the correlation function by a nonzero spectral 
width, section 5 contains an example with semi-realistic pion spectral function in 
hot nuclear matter. Finally we draw some conclusions and discuss the experimental 
relevance. 



2. Outline of the method 

In "ordinary" quantum mechanics, a statistical state of a quantum system is de- 
scribed by a statistical operator (or density matrix) W, and the measurement of 
an observable will yield the average 

\ TT[£it,x) W] 
^(''^7^ Tr[W] ' 

where the trace is taken over the Hilbert space of the quantum system and £ is 
the hermitean operator associated with the observable. In thermo field dynamics 
(TFD), the calculation of this trace is simplified to the calculation of a matrix 
element 

with "left" and "right" statistical state defined in terms of the two different com- 
muting representations (see refs. |l^ for details). 

In this state, we consider a complex, scalar boson field describing spinless 
charged excitations in a statistical system not too far from equilibrium. In the 
spirit of the first remark, one could think of this field as describing positive and 
negative pions in nuclear matter. 

According to the reasoning above, this "thermal" boson field is described by 
two field operators (f>x, 4>x and their adjoints (/)*, (j)*, with canonical commutation 
relations 

[<f>{t,x),dt(f>*{t,x')] = iS^{x-x') 
^{t,x),dt^*{t,x')\ =- iS^ix-x') (3) 



but commuting with each other. These two fields may be combined in a statistical 
doublet, see r ef. |l| for details. 



The free as well as the interacting scalar field can be expanded into momentum 
eigenmodes 

^^'1= (4_(t)e~"=- + a,+ We'^-) 



(4) 



iEt 



k is the three-momentum of the modes, therefore in this notation a\_{t) creates a 
negatively charged excitation with momentum fc, while afc+(i) annihilates a positive 
charge. Henceforth the two different charges are distinguished by an additional 
index I = ± whenever possible. 

For the free case the commutation relations of the a-operators at different times 
are simple, while they are unknown for the interacting fields. However, we want to 
go only one step beyond the free field case, i.e., we approximate the fully interacting 
quantum fields by generalized free fields ]lO[ |. In this formulation, the operators a,a 
do not excite stable on-shell pions. Rather, they are obtained as an integral over 
more general operators ^, ^ with a continuous energy parameter E [T5| p^ : 

) = ldEjd^qA]'\E^k) {BT\E,q,k)y ( ) e~ 

-Hjn"^ = IdE ld^qA\l\E,k){ 1*^' ^ 6Ki^,g,fe)e'^', (5) 

afcK^j / \ "KEql j 

where ,B is a 2 x 2 matrix, the weight functions AiiE^k) are positive and have 
support only for positive energies, their normalization is 

^dE E Ai {E, fc) = i ^dE Ai {E, k) = Zu . (6) 



The principles of this expansion have been discussed in rcf. ||ic|] , its generalization 
to non-equilibrium states was introduced in ref. ||l^. For equilibrium states the 
combination 

A{E,k) ^ A+{E,k)Q{E) - A-{-E,-k)Q{-E) (7) 

is the spectral function of the field ipx and the limit of free particles with mass m 
is recovered when 

A{E, k) — > sign(£;) - k^ - m^) = sign(i;) S{E^ ~ ujI) . (8) 

For non-equilibrium systems, the existence of a spectral decomposition cannot be 
guaranteed [p^ . We may expect however, that close to equilibrium the field prop- 
erties do not change very much. Thus, with this formalism we study a quantum 



system under the influence of small gradients in the temperature, with local spec- 
tral function A{E, k). Corrections to such a picture only occur in second order of 
temperature gradients jl^, |l^ . 

A thorough discussion of the 2x2 Bogoliubov matrices was carried out in ref. 
[12|. For the purpose of the present paper, we simply state their explicit form as 



Bi{E,q,k) 



{S^q-k)+Ni{E,q,k)) -Ni{E,q,k) 
-5^{q-k) 5^{q-k) 



(9) 



where N{E, q, k) is the Fourier transform of a space-local Bose-Einstein distribution 
function 



Ni{E,q,k) 
ni{E,z) 



1 



(27r)3 



jd^ze-'^t-^^^ni^E.z) 



1 



g/3(2)(£;-Mi(2)) _ 1 



(10) 



Here we have assumed a distribution function that only depends on the energy pa- 
rameter E and on the space coordinate z. The expansion allows for a generalization 
of this, to more general distribution functions depending also on the momentum 
{q + k)/2. 

We have argued, that our ansatz for the fields gives rise to a local spectral 
function. A moving particle however feels an influence also of the gradients of this 
local equilibrium distribution. Consequently also the propagator for the fields we 
consider is correct beyond a local equilibrium situation, to be precise it is correct 
to first order in the gradients of ni{E, z). 

To complete the brief description of the TFD formalism, we specify the com- 
mutation relation of the various operators in our expressions. The ^-operators have 
commutation relations 



^5u'5{E~E')5^{k~k') 



(11) 



Similar relations hold for the ^ operators, all other commutators vanish, see pO[ |. 
It follows from these definitions, that 



aki{t),a\,i,{t) 



Zki 5u' 5^{k-k') 



aki{t),al,i,{t)\ = Zki5u'6\k-k') 



(12) 



are the equal-time commutation relations for the a, a operators. 

The ^, ^ operators act on the "left" and "right" statistical state according to 

iEM\\W)) = Q, lEki\\W))=Q, ((l||eL=0, mltki^Q yE,kJ = ±l. (13) 



With these rules, all bilinear expectation values can be calculated exactly. Higher 
correlation functions have a perturbative expansion in the spectral function. 



3. Two-particle correlation function 



Of the higher eorrelation functions, we are interested in the two-particle correlation 
function, which is the probabihty to find in the system a pair of pions with mo- 
menta p and q. For the non-equihbrium system we are considering, this correlation 
function is 

{ali(t)al,{t)a,v{t)a,i{t)) T{p.,q)Hq,p) ,,,, 

For simplicity, we abbreviate the mean momentum of this pair by Q = (q +p)/2. 
The function !F{q,p) is calculated using the standard rules of thermo field dynamics 
given above. One obtains 

oo 

:F{p,q) - JdE jd?z {Ai{E,p)ME,q)y e'^P-"^^ ni{E, z) , (15) 



where one may also insert a ^-dependent spectral function without violating the 
accuracy to first order in the gradients. 

How these gradients enter the above expressions may be seen when performing 
an expansion of A around the mean momentum Q. 

^ip,q) = ^°(p,g) 

oo 

+ JdE J d^z Sp-i> {WqAi {E, Q, z) V^ni {E, z)) 



+ 0{^ln). (16) 
Here, the lowest order term 

oo 

^{P.q)^ JdE Jd''zc'^''-P'>^Ai{E,Q,z)ni{E,z) (17) 



is the local equilibrium contribution to the J- we have obtained above. 

For a possible generalization, i.e., to explicitly momentum dependent ni, it is 
worthwhile to note that the gradient term in (^6|) is just one half of the Poisson 
bracket of A and n [|l^. Furthermore, we find that T{p,p) = J-''^{p,p), i.e., the 
denominator of the correlation function is not affected by the gradient expansion. 

We therefore obtain as the local equilibrium two-particle correlation function 
the expression 

Joe. , . .F°(p, q)T'^{q,p) 



However, the full non- equilibrium correlation function cw (p, q) is the one measured 
experimentally. 

The exact correspondence between the local equilibrium result and other cal- 
culations of the correlator ^, ^, Q may be found when inserting the free spectral 
function from (H): 



J d^zni{ujp,z)j (^J d^z' ni{ujq,z') 



(19) 



Obviously, one may not insert the free spectral function into equation (^^ for J^. 
This is not a flaw of the derivation, but suggests - as expected - that the limit of 
zero spectral width at finite temperature is ill-defined [p^. 



4. Simple spectral function 

In this section we study the difference between the local equilibrium correlation 
function (|l^) and the non-equilibrium result ( |l5| ) in more detail. To this end we 
calculate the correlation functions with a simple parameterization of a boson (pion) 
spectral function, 

-^'(^'^)-?(^^3T^fFTlW ^''^ 

where flp — ^ rn^ + + 7^ and m-,^ — 140 MeV. To gain information about the 
maximal influence exerted by the occurrence of a nonzero spectral width, we use 
an energy and momentum independent 7 equal for both charges. 

The temperature distribution is taken as a radially symmetric gaussian, 

T(^)=r(r)=Toexp(^-^^ , (21) 

with chemical potential = and Rq — 5 fm. 

The local equilibrium pion distribution for a given momentum k is obtained by 
folding z) with the spectral function. Hence, the mean radius of the particle 
distribution function acquires a 7-dependence. We define the rms radius orthogonal 
to the direction of k as 



Rr 



J- I J = Jdrr^ JdEA{E,k) (e^/^M - l) \ (22) 







Note, that i?rms is not the 3-dimensional rms radius of the distribution function 
(which would be 14/12)- Rather, Rrms is half the product of angular diameter and 



distance between detector and source. A constant temperature over a sphere of 
radius Ro would yield an i?rms = Ro/VS, while its 3-D rms radius is i?o-\/3/5. 

In fig. |l] we have plotted the correlation functions for two different constant 
values of the parameter 7. Clearly, for 7=50 MeV the correlation function cw {p, q) 
agrees quite well with the local equilibrium result, and very closely resembles a 
gaussian. 

However, for smaller 7=5 MeV the non-equilibrium correlation function is much 
narrower in momentum space than the local equilibrium result, it also deviates from 
a gaussian form. Nevertheless we may approximate it by such a simple functional 
form in order to extract quantitative information, i.e., 

cw{p,q)^l + exp{-R\p-q)^) (23) 

and similarly for c}°*^(p, q) with parameter R[q(.. 

In figure H we show the two fit parameters R and i?ioc function of 7, together 
with the 7-dependent rms radius of the particle distribution. We find that the 
measured correlation radius R is always larger than Riq^. , with a minimum reached 
at 7 « 26.7 MeV. 

The plot may be divided in two regions, with a boundary at 27i?rms[7] = 1 <^ 
7 « 30.2 MeV. For these two regions we find 

27i?rms < 1 i? ~ i?ioc + > ^rms > i?ioc ^24) 

27-Rrms ^1 ^ R~ R\oc Rrms > Rloc ■ 

For larger 7, the small differences between R, Riqq and i?rms may be attributed 
to our use of a gaussian temperature distribution: n[T{r)) is not strictly gaussian, 
only in the (unphysical) limit 7 — > 00 one reaches R — R\oc = Rrms — Ro/^f^ 

The interpretation of this result is straightforward: A finite lifetime or nonzero 
spectral width 7 > of the bosons inside the source is essential, if one wants to 
infer the thermal source radius Rq from correlation measurements. To be more 
precise, only for 27i?rnis > 1 the correlation function measures the mean diameter 
of the particle distribution function n{E,z). 

This result is in agreement with our view of the relaxation process of a non- 
equilibrium distribution function: The relaxation rate is, to lowest order, given by 
the spectral width of the particle . Consequently, zero 7 corresponds to a system 
without dissipation. In such a system, only quantum coherence effects exist, and 
thus the correlation function approaches the "quantum limit" 

\im cw{p,q) ^ 1 + SwSpq . (25) 

For this case, the correlation radius obtained by a gaussian fit becomes infinite. 
We may also view, for a given energy, I/7 as a measure for the spatial size of the 
pion "wave packet" , which must be smaller than the object to be resolved. In other 
words, the mean free path of the bosons must not exceed the object size to produce 
correlations. 
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0.14 GeV 0.938 GeV 


1.232 GeV 


0.12 GeV 


0.155 fm 



Table 1: Coupling constants and masses used in the calculations of this work. 
The value of g' was chosen to allow for a direct comparison with simulation codes for heavy- 
ion collisions, a more realistic value to describe pion scattering data would be g' « 1/3. 



5. Semi-realistic spectral function 

To get a more realistic result for the non-equilibrium two-pion correlation function 
measured in the thermal radiation from a hot spot in nuclear matter, we use the 
the spectral function derived in ||l^, |l^. It includes the coupling of pions to Ass- 
resonances in nuclear matter, which are taken to have a constant spectral width T 
by themselves. 

It was argued in ref. p2[ , that using such an approximate spectral width for 
the A33 resonances constitutes the only way to achieve an analytical solution of 
the A-hole polarization problem in nuclear matter. An even more realistic energy- 
momentum dependent spectral width for the A33 resonance can be treated only in 
a fully self-consistent numerical treatment involving dispersion integral techniques. 

To first order in such a constant F, the pionic spectral function is 

^ _ '^^'^(-^) TEuoA 

Note, that this function is coordinate dependent. The energies in the denominator 
are 

^± = \ [eIa + (r/2)^ +El± ^{El^ + (r/2)2 - Elf + Ak'C{z)u:^ , (27) 
with functions 



UJA = EA{k) - Mn = Jk^ + Ml- AfjN 



O / XTT \ / 1 

^ / Jna \ / ^0 ^ ^0 



ENAik) = V^A(fc) {iVAik) + g' Ciz)) . (28) 

The baryon number in each small volume and hence the baryon density is a constant 
parameter of the calculations, for free particles with bare "on-shell" energies 



ENip) = + ^^N and EAip) = Jp^ + Ml (29) 



this baryon density is obtained as 

= ^ J-^^NiENip),z) + 16 J-^nA{EA{p),z). (30) 
The distribution functions are taken as local Fermi-Dirac functions 

^) = e(E-,..^)/Ti.) + 1 ' (31) 

with temperature T{z). As temperature distribution we use the same as in the 
previous section, eq. (|2l|), but with different central temperatures Tq. 

Since the temperature depends on the spatial coordinate, fixed baryon density 
implies that the "baryochemical" composition of the hot spot changes with coordi- 
nate z. In the center of the hot spot baryons are, to a large extent, present in the 
form of A-resonances. Outside the hot region baryons are "only" nucleons - and 
as shown before, the pion properties there do not influence our calculation. 

In fig. 1^, we have plotted the correlation function for a pair momentum of 
Q = (p + q)/2 = 100 MeV at the two values To= 100 MeV and To = 160 MeV. In 
this momentum region, the spectral width of the pion is small due to its pseudo- 
vector coupling with baryons. This we assume to be a general feature of pions 
in nuclear matter, although one may argue about the exact value of 7. We find, 
that at higher temperature the non-equilibrium correlation function (which we had 
assumed to be the one measured experimentally) is much narrower in momentum 
space than the local equilibrium function. 

The local equilibrium correlation function however is close to the the rms- 
radius of the thermal source. The correlation radii obtained by a gaussian fit to 
the non-equilibrium as well as the local equilibrium distribution are given in table 
^ Following these results we conclude, that the effect we propose is absent at 
higher momentum of the pions, where the p-wave coupling to nuclear matter is big 
enough to give it a sufficiently large spectral width for local equilibration. In the 
low momentum region, the measured correlation radius overestimates the source 
size by as much as 30 - 40 %. 



T (p + g)/2 R Rrms R/Rrms 



100 MeV 


100 MeV 


3.99 fm 


2.89 fm 


3.14 fm 


1.27 


160 MeV 


100 MeV 


4.86 fm 


3.23 fm 


3.51 fm 


1.38 


100 MeV 


350 MeV 


3.23 fm 


3.16 fm 


3.34 fm 


0.97 


160 MeV 


350 MeV 


3.57 fm 


3.49 fm 


3.73 fm 


0.96 



Table 2: Correlation radii obtained by gaussian fit to the correlation function with 
semi-realistic spectral function. 



6. Conclusions 



Before we draw a final conclusion from our work, we must emphasize that it is 
still too early to use our result for the correlation function in a direct comparison 
to experimental data. For any realistic situation, we certainly have to take into 
account also the partially coherent production of pions: A substantial fraction of 
pions arriving in a detector stems from the free-space decay of A33 resonances, thus 
forcing cu'{p,p) < 2. 

Furthermore it seems worthwhile to note that our derivation does not contradict 
existing work on the correlation function. In an early semiclassical treatment, an 
infinite lifetime (equivalent to 7 = 0) of an excitation was found to produce a 
correlation function c = whereas we find an infinitely narrow peak in this 

"quantum limit" , see eq. ( |25| . 

Another paper based on the Schwinger-Keldysh method also finds that the 
total production rate of particles is proportional to an energy integral over the off- 
diagonal self-energy components, i.e., to the integral over 7 • n in our formalism 

However, due to their explicit quasi-particle approximation the authors of ref. 
[7\ lose the consistent treatment of the non-equilibrium effects to first order in the 
temperature gradients. 

Let us briefly remark on the physical situation present in stars emitting pho- 
tons: For these we may assume a thermal width of the photon which is very small, 
e.g., 7 « aT « 0.5/137 eV. However, the source radius i?rms of a star is, in general, 
so big that the condition 27i?rms ^ 1 is always satisfied. Hence our results do not 
affect correlation radii measured for astronomical objects. 

Having in-lined our calculation with existing work, we may now carefully con- 
clude the following physical effect relevant for relativistic heavy-ion collisions: When 
pion pairs with a relatively low momentum are created in hot nuclear matter, their 
p-wave interaction with the surrounding medium is small. Hence, in their move- 
ment outwards from the hot zone they do not have a sufficiently short mean free 
path to be locally equilibrated. Thus, their correlation length is dominated by their 
mean free path (w I/7) rather than by the thermal source radius. 

As we have shown, this leads to a correlation function which is narrower in 
momentum space than expected in a local equilibrium situation: Compare the solid 
and the dashed curves in fig. ^. Consequently, measuring the correlation function 
of low momentum pions in a non-equilibrium state overestimates the source radius. 
For a semi-realistic pion spectral function and a pair momentum of 100 MeV, we 
find this effect to be as large as 30 - 40%, depending strongly on the actual source 
size. 

Turning to experimental results of NA44, we learn that correlation measure- 
ments of pions and K-mesons emanating from relativistic heavy-ion collisions yield 
comparable fitted correlation radii of 3-4 fm ||] . As we have shown, it might be pre- 
mature to conclude from these measurements that the source size for both mesons 
is similar: One would have to compare the mean free path of both particle species 



before such a conclusion. 

More recently the experiment NA49 has measured pion correlation functions 
in central Pb+Pb collisions at 33 TeV. Resulting was a correlation radius of 6-7 
fm, as compared to 4.5 fm in S+Au collisions ||^. It is a particularity of our semi- 
realistic pion spectral function that it is narrower for higher temperature and low 
momenta. Physically, the first effect stems from the higher A33 abundance in the 
hotter system: A-particle/nucleon-hole pairs are less easily polarized. The second 
effect is due to the p-wave coupling of pions to nucleons, as was pointed out above. 

Consequently, the effective 7 of the pions drops with increasing temperature - 
and one measures a higher correlation radius. If we presume this effect, which was 
deduced for equilibrated matter, also to hold in the highly non-equilibrium situation 
present in the experiments, we can state that the higher correlation radius of the 
NA49 data might indicate a higher temperature of the reaction zone rather than 
its bigger size. However, this interpretation depends crucially on the momentum 
range where the pions are measured. The uncertainties in this statement indicate 
the importance of calculations for mesonic spectral functions in hot nuclear matter. 

The effect of narrowing the correlation function might be even more pronounced 
due to the short time-scales of relativistic heavy-ion collisions, which may suppress 
the spectral width (~ collision rate) of pions even more. Indications for such a 
behavior are obtained in simulation codes, where at high enough collision energies 
one may reproduce experimental data with free particle collision rates. Thus, in- 
stead of speaking about "temperature" , which always is connected with the notion 
of a partial equilibrium, we may also reformulate our conclusion even more sharply: 
Measured correlation functions become narrower due to non-equilibrium effects in 
statistically radiating sources. A bigger correlation radius therefore is an indica- 
tion, that a system of colliding heavy-ion collisions is farther from local equilibrium 
in the collision zone. 

Finally we point out that in common calculations of the correlation function 
one has to introduce ad-hoc random phases between several classical sources and 
then obtains only the local equilibrium correlator c}j^'^(p, q). Instead, we relied 
on a proper field theoretical treatment which incorporates non- equilibrium effects 
correctly up to first order in gradients of the temperature. We found, that the 
non-equilibrium character of the system must be taken serious when calculating 
the correlation function - which answers the secondary question we asked in the 
introduction. 

Coda 

A large part of the work of E.Wigner was dedicated to the study of symmetries in 
physical systems. The fact, that particles in a thermal state must have a nontrivial 
spectral function rather than a sharp " mass-shell constraint" , is connected to the 
breaking of such a symmetry: A thermal state has a preferred rest frame and 
therefore violates the Lorentz invariance of the usual field theoretical ground state 
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Figure 1: Correlation function for a pion "hot spot" with temperature Tq^IOO 
MeV. 

Simple spectral function from eq. (po[), (p + q)/2= 100 MeV. 
Thin lines: 7=50 MeV, Thick lines; 7=5 MeV. 
Solid lines: non-equilibrium correlation function 
Dashed lines: local equilibrium correlation function ([L8[). 
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Figure 2: Correlation radius of a pion "hot spot" with temperature Tq; 
Simple spectral function from eq. @, {p + q)/2= 100 MeV. 
Thick solid line: gaussian fit radius R of the correlation function ([l4|). 
Thick dashed line: gaussian fit radius -Rjqq of the correlation function (p^. 
Thin solid line: 7-dependent rms radius -Rrms- 
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Figure 3: Correlation function for pions from a "hot spot" with different central 
temperature. 

Semi-realistic spectral function from eq. (p^), (p + q)/2= 100 MeV. 
Top panel: Central temperature Tq — 100 MeV, 
Bottom panel: Central temperature To = 160 MeV, 

SoUd thick line: non-equilibrium correlation function (wj). 

Dashed line: local equilibrium correlation function (^). 

SoUd thin line: gaussian with correlation radius equal to -Rrmo (see table). 



